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ABSTRACT 


Observations show that hard X-ray burst and UV lines rise and fall 
simultaneously on time scales of seconds. Hydrodynamic simulations of beam- 
heated atmospheres, based on collisional transport, however, produce only a 
gradual fall in UV emission, when the beam flux falls, due to the long time 
scale of conductive relaxation. It is suggested that this discrepancy might 
be explained by onset of plasma turbulence driven by the strong heat flux or 
by the beam return current going unstable. Such turbulence greatly reduces 
electrical (a) and thermal (k) conductivities. Fall in a reduces the hard 
X-ray flux by enhanced ohmic dissipation of the return current, while fall 
in k may cause the UV line to fall by reducing the transition region thick- 
ness. 


1 Introduction 

Simultaneous SMM data in hard X-ray and UV lines exhibit synchronism of 
both rise and fall of impulsive spikes down to time scales of order of seconds. 
The driving mechanism for both is believed to be a thick target electron beam. 
However, numerical simulation of the atmospheric response to such a beam 
(Mariska and Poland, 1985) does not reproduce synchronous fall in the two 
emissions. In such simulations, UV line emission rises because of the direct 
heating effect of the beam and because of the increased conductive flux in the 
transition region driven by beam heating of the corona (cf Section 2) . Though 
the direct heating turns off with the beam, the fall in the enhancement due to 
conduction is gradual because the corona acts as a heat reservoir. 

This paper gives a very preliminary discussion of how onset of beam 
generated plasma turbulence might modify the simulation results through 
anomalous transport effects. 
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2 Line Formation in the Flare Transition Region 


An injected electron beam dominates flare atmosphere heating in the 
corona and in the low chromosphere, where the temperature gradient is rather 
small, while in the upper chromosphere and transition region, after a brief 
initial transient, the steep temperature gradient results in energy transport 
being dominated by thermal conduction (driven by the coronal beam heating) - 
e.g. Shmeleva and Syrovatskii (1973) and subsequent authors. In a quasi- 
steady transition region there is negligible direct beam heating, because of 
the small column mass, and essentially constant conductive flux down to the 
temperature (2l 10 5 K) where radiative losses become stable and maximal, and 
where densities become high enough for conductive deposition to be radiated 
away. Furthermore, because of the small column mass, the transition region 
can be adequately described as having spatially uniform pressure. 

For a conductive energy flux F emanating from the corona, the thermal 
structure of the transition region is then described by 

k(T) ~ = F (1) 

where k(T) is the thermal conductivity at temperature T, at geometric depth z. 
Since the bulk of the energy of a beam is deposited in the corona, F is 
essentially the energy flux of the beam. 

An optically thin collisionally excited transition region line of wave- 
length X, with emissivity n 2 f^(T) (erg cm“ 3 s -1 ) at plasma density n, 

temperature T, will have a total luminosity per unit area 


I 


X 


n 2 (z) (T (z) ) dz 

Z 


( 2 ) 


where n is determined by the constant pressure condition nT = 


n o T o 


where n r 


is the atmospheric density at the upper chromospheric level where the bolo- 
metric radiative loss curve maximises at temperature T Q Ck 10 5 K) . n Q will, 
of course, depend on the magnitude of F, the higher value associated with a 
flare pushing the base of the transition region deeper in the atmosphere to 
the depth where n Q is high enough to radiate off the input. For lines formed 
in the upper transition region (T > 10 5 K), (2) then becomes 
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As a first approximation, a typical line may be approximated as being formed 
over a fixed interval AT centred on the peak temperature T of line formation 
(3) can be approximated, using (1), as 


and 


*x = 


f X(T) dT 


T 2 F/k (T) 
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where k(T) is a mean value of k over AT approximated by its value at line 
peak T. 

It follows that if k(T) remains of the same form, 1^ is determined by 
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n 2 /F. Increase in F alone during a flare would result in a decrease in I, 
o 1 11 1 A 

due to the transition region becoming thinner. However, as already noted, 
increase in F also results in increase in n Q and steady state numerical 
simulations (Emslie, 1985) show that the rise in n Q more than offsets the 
decrease in AZ = k(T)AT/F and the nett result is a rise in UV flux 
accompanying an increase in thick target hard X-ray flux due to the increased 
F. Time dependent simulations (Mariska and Poland, 1985) confirm this and 
show that the rise is closely synchronous but that when F falls, the thermal 
conduction relaxation time is too long for 1^ to track F. All of these 
calculations, however, assume that k(T) retains the same form throughout and 
in particular neglect the consequences of onset of anomalous transport 
processes. 


3 Effect of Onset of Plasma Turbulence on Transition Region Lines 
3.1 Qualitative Description 

Suppose that instead of assuming the beam flux F to peak at the observed 
peak of a radiation spike, while transport is still classical, we suppose 
that during its rise F exceeds a threshold for anomalous transport effects to 
set in. This may be due to one or more of three processes: two-stream 
instability of the beam itself, resulting in Langmuir wave generation (e.g. 
Emslie and Smith 1983; McClements et al 1985 - this Workshop); drift current 
instability of the beam driven return current resulting in generation of 
electrostatic (ion-cyclotron, ion-acoustic) waves (e.g. Hoyng et al 1977; 
Cromwell ejt al 1985, Holman 1985 - this Workshop); heat flux instability 
resulting in electrostatic wave generation by the currents associated with 
steep thermal gradients (Mannheimer 1977, Brown et_ al_ 1979, Smith and 
Lilliequist 1979). Here we will concentrate our attention on the last two 
processes as it is these which will directly affect thermal and electrical 
conductivities . 

When such plasma waves appear, the effective electron collision 
frequency will increase, resulting in a fall in both electrical (a) and 
thermal (k) conductivities. Decrease in a will increase the electric field 
required to drive the return current, reducing beam electron lifetimes, and 
consequently the X-ray bremsstrahlung from the beam. Decrease in k, on the 
other hand, will result in a decrease in the transition region thickness (1) 
and consequently in 1^ provided n Q does not rise enough to offset this factor. 
The extent to which n Q will eventually rise for a given change in k will 
have to be determined by a steady state numerical calculation (cf Machado and 
Emslie 1979). For the moment we will assume that n Q is determined mainly by 
the input flux F rather than by details of k(T) and that the k factor in 1^ 
overwhelms the n 2 factor. Likewise, the time scale on which the transition 
region thins wil? have to be found by time dependent numerical simulation. 

We note, however, that a change in structure purely due to a change in 
coronal heating when F changes will take roughly the time of propagation of 
a thermal front along the coronal loop length. On the other hand, re- 
adjustment of the transition zone structure due to an in situ change of k 
( typically in a few plasma periods) will occur in the much shorter time needed 
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for a thermal front to cross the transition region. 


3.2 Conditions for Onset of Anomalous Transport in the Flare Transition 
Region 

If beam electrons are injected into the transition region at a total 
rate ^(s -1 ) over an area A then the condition for return current instability 
can be written 

n (^) /2 f(T e /T i> < ^ /A (5) 

' e ' 


where f declines from about unity as T e /T. increases from unity (e.g. Hoyng 
e£ al_ 1977, Brown and Hayward 1981). 

On the other hand, the heat flux F = E &/A is likely to drive wave 
generation if it significantly exceeds the saturated value ^ -^r / kT e Y ^t or 
roughly “e ' 

(kT e \V2 kT e 

n ra — £ -? /A < 6 > 


where E is the mean energy of beam electrons. 


From (5) and (6) we see that which effect sets in first depends on the 

value of f compared to kT„/E . With E > lOKeV and kT < 6.5KeV for the 

^ ** 0 ~ 

transition region, it seems almost certain that heat flux instability will 
set in first except for very high T £ /T£ values and near the top of the 
transition zone. (Noting the condition nT £ = constant we see that the left 
side of (5) varies as T'V 2 f (TjT/) which increases with depth in the 
atmosphere whereas the left side of (6) varies as T # and so decreases with 
depth - i.e. the heat flux is most unstable at low T .) In absolute 
numerical terms, the separate criteria are as follows. 


Heat Flux Instability 

If the transition region pressure is P(dyne cm -2 ) = 10 2 P ? % nkT 
with T = 10 6 T 6 (K), >= 10 36 ^ 36 (s _1 ), A = 10 18 A 18 (cm 2 ) * 

then (6) becomes 


and 

and "e = 10E 1 (KeV) 


g36El > 2 3P T V2 
A 18 ' 26 


(7). 


Thus for typical flare transition region pressures P 2 24.0.3-3 (e.g. 

Machado et al 1979) and temperatures T & 1 and with the electron beam 
parameters ^ ^ E x ^ 1 typical of large hard X-ray bursts, heat flux 

instability in the transition zone is likely if the injection area 
A <_ 4 x lCr'cm 2 which is very possible. 
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Return Current Instability 

With the same notation and f ■ 0.1 f_^, (5) becomes 

9 . 


IT- 1 2 'V-i 
a 18 


( 8 ). 


This condition will be satisfied only if f is particularly small, due to 
large T e /T^, or if the injection area A is << 10 18 cm 2 , for typical flare 
transition region pressures. It is more readily satisfied in the corona, 
where P is smaller (and larger) , and so may contribute strongly there to 
beam deceleration, and hence to the fall in bremsstrahlung from the beam. 


3.3 Quantitative Effect of Plasma Turbulence in the Transition Region 

Both the electrical (a) and thermal (k) conductivities are determined 
by the effective collision frequency v e ^ of the electrons, viz 

< 6 x lO^n^T^/v ^ ergcm -1 s^K -1 (9) 

and 

a 2 x 10 18 n- /v 


.-1 


10' eff 


( 10 ) 


where n 


n 


10 10 10 ‘ 

In the classical regime v e f£ ■ V coll §* ven by 


v coll ■ 7 ’ , l° a "10 /T 6 % 


(ID 


whereas in the presence of (e.g. ion-acoustic) waves of energy density W, 

( 12 ) 


v rr becomes of order 
eff 


eff 


V (ds) ■ lo No w (=) 


where v^ e is the plasma frequency. Thus the classical and anomalous transport 
coefficients become 


and 


K . “ 10 9 T, % 

class 6 

(13a) 

"an ■ 6 * / (sf) 

(13b) 

a , = 3 x 10 15 T,^ 

class 6 

(14a) 

“am - 2 “ “No 14 / (Jr) 

(14b) 


If then, during a UV/HXR burst, waves are generated by either return 
current or heat flux instability, at wave onset the thermal conductivity and 
hence the transition region scale thickness will fall by a factor (Equation 
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(l)) 


AZ 

"az 


'AN 

class 


AN 

C class 


* 6 * 10 ' 7 “io ft V %/ (55) < 15 > 


The actual value of W/nKT depends on how the waves are driven. Simulations 
of return current instability (Cromwell e£ al 1985 - these proceedings) lead 
to values of W/nKT as high as 10 ^ at saturation and around 10~^ in marginal 
stability. In the more relevant case of heat flux instability. Brown et al 
(1979) found that in the case of extreme temperature gradients, unstable 
generation of ion sound waves led to W/nKT 2 l (m /m )^ 1/43 for which (15) 

gives e p 


AK 


AK 


AN 

class 


i 2 x 10 


i-4 p 


/t 6 2 


Consequently, neglecting any compensating factor due to increase of P Q , 
Equations (4) and (16) show that a small increase in F past a critical thres- 
hold can result in a reduction in UV line strength by several orders of 
magnitude. Furthermore, we expect this decrease to be very fast - roughly 
the time it takes for F to rise above the instability threshold across the 
temperature domain of formation of the line concerned. 


Whether wave generation is driven by heat flux or return current instabil- 
ity > its effect on a will act on the electron beam in two ways. Firstly, 
it will increase the electric field driving the return current and so 
decelerate the beam in the region of its propagation. Secondly, and more 
speculatively, it may feed back on the region of acceleration of the beam and 
interfere with beam production (Brown and Melrose 1977). To see the 
consequences of the first effect on the HXRB, we can compare the thick target 
bremsstrahlung yield in the case of collisional losses only (the usual thick 


target case - Brown 1971) with that when strong return current losses 
dominate. In each case the bremsstrahlung photon yield at energy e from an 
electron of initial energy E is roughly v = nQ^vx where is the 

bremsstrahlung cross section (appropriately averaged over electron energies 
between e and E ) , v is the electron velocity and x is stopping lifetime. 
In the classical (collisional case) x = /Zire^Anv and so 

Brems Qb^c* 2 

— nr- 


v , 'V 
class — 


2ire 4 A 


(photons per electron) 


(17). 


In the case of return current losses dominating (e.g. Brown and Hayward 1981) 
T = m e v o /(e j /a > = oAE Q /e 2 Jv o so that the photon yield is 

Brems 


AN 


^ n Q B a AN AE o /e2 ^ 


(18). 


Taking the ratio of (18) to (17) substituting, for from (14b) and insert- 

ing numerical values then gives 

'v. - - % 




i 3 x io 


-8 


Brems 


n 10 A 18 / / W \ 

^ 36 E 1 ' ^ nKT ^ 


(19) 
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which even with W/nKT as small as 10** 1 * implies a fall in thick target 
bremsstrahlung yield of over 3 orders of magnitude, effectively switching 
off the hard X-ray production entirely. 


4 Conclusions 

In summary, when the electron beam flux in a beam heated flare becomes 
large enough for onset of wave generation either directly by the beam return 
current, or via production of heat flux beyond saturation, we expect an 
immediate reduction in both hard X— ray burst intensity (due to anomalous 
return current dissipation) and simultaneously in the emission of UV lines 
from the transition region (due to reduction in its thickness by anomalous 
thermal conductivity). The exact magnitude and time scales of this effect 
require further investigation by numerical simulation, hopefully as a sequel 
to this Workshop. In particular, the effect of reducing k on the transition 
region pressure (i.e. n Q ) and its effect on the above conclusions is the topic 
of Workshop related collaborations. 
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